ApJ, 486, in press 



On the Origin of the UV Upturn in Elliptical Galaxies. I. Sensitivity of UV 
Population Synthesis to Various Input Parameters 

Sukyoung Yi 1 

Department of Astronomy, Yale University, P.O. Box 208101, New Haven, CT 06520-8101, 
NASA/Goddard Space Flight Center, Code 681, Greenbelt, MD 20771 
email: yi@shemesh.gsfc.nasa.gov 

Pierre Demarque and Augustus Oemler, Jr. 2 
Department of Astronomy, Yale University, P.O. Box 208101, New Haven, CT 06520-8101 

demarque@astro.yale.edu oemler@ociw.edu 

ABSTRACT 

We present models of the late stages of stellar evolution intended to explain the 
"UV upturn" phenomenon in elliptical galaxies. Such models are sensitive to values of 
a number of poorly-constrained physical parameters, including metallicity, age, stellar 
mass loss, helium enrichment, and the distribution of stars on the zero age horizontal 
branch (HB). We explore the sensitivity of the results to values of these parameters, 
and reach the following conclusions. 

Old, metal rich galaxies, such as giant ellipticals, naturally develop a UV upturn 
within a reasonable time scale - less than a Hubble time - without the presence of 
young stars. The most likely stars to dominate the UV flux of such populations are 
low mass, core helium burning (HB and evolved HB) stars. Metal-poor populations 
produce a higher ratio of UV-to-V flux, due to opacity effects, but only metal-rich stars 
develop a UV upturn, in which the flux increases towards shorter UV wavelengths. 

Model color-magnitude diagrams and corresponding integrated spectra (for various 
values of age, metallicity, helium enrichment, mass loss efficiency, initial mass function, 
and the HB mass dispersion factor) are available on S.Y.'s world wide web site 



http: / / shemesh.gsfc.nasa.gov/model.html , 
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1. Introduction 



The evolutionary population synthesis (EPS) technique, first introduced by Tinsley (1968; see 
[Tinsley 1980 and Bruzual & Chariot 1993 for review), has its fundamental basis in stellar evolution 
theory. Unlike the static population synthesis^ technique (e.g. McClure fc van den Bergh 196% 
Spinrad fc Taylor 1971fc |Faber 1972j ; IQ'Connell 1976] ; pickles 1985b] ; see |Q' Conner! 1987] for review) 



that adopts physical constraints from the observational data, EPS relies on theoretical constraints, 
most importantly the evolutionary time scale of stars. 

Some significant advantages of EPS include the following: (1) the theoretical results can 
provide physical understanding when the results fit the observations reasonably because the 
solution is unique, (2) it allows us to understand not only the present observable quantities but 
also the evolution of those quantities, (3) it is easy to implement the theoretical predictions into 
the code, and (4) it is easier to investigate the sensitivity of the result to input assumptions that 
are used in the modeling. While aspects (1) and (2) have been widely recognized, the significance 
of aspects (3) and (4) is still poorly understood. 

As an example of aspect (3), EPS enables us to include hypothetical stars that are predicted 
to exist and significant contributors to the integrated flux but rarely observed. For instance, 
several stellar evolution theory groups have predicted that, under certain conditions, low mass 
core helium-burning stars (horizontal branch - HB - stars) become UV-bright before they become 
white dwarf (WD) stars ( pemarque fc Pinsonneault 1988] ; |Grcggio fc Renzini 199Cl| ; Castellani fc 



Tornambe 1991| ; Horch, Demarque, & Pinsonneault 1992). These UV-bright stars, the so-called 
slow blue phase (SBP) stars ( Horch et al. 1992[ ) which include both AGB-manque stars ( Grcggio] 
fc Renzini 1990] ) and post-early-AGB (PEAGB, [Castellani &: Tornambe 1991| ) stars, can be 
realistically taken into account only through the EPS technique. Yi, Demarque, & Kim (1997, 
hereafter YDK) provided a detailed physical description of the UV-bright phase. 

On the other hand, the reliability of EPS is restricted by uncertainty in the stellar physics. 
The SBP stars give an example again. They were unknown until recently. Thus, early EPS studies 
were not able to include them, despite their probable significance in the UV study. However, 
such pure theoretical predictions from stellar astrophysics may not be realistic: such stars may 
not exist. Then, the EPS using such theoretical constraints would result in wrong conclusions. 
Therefore, in earlier times when many of the stellar evolutionary phenomena other than the 
main sequence (MS) and red giant branch (RGB) were unclear, the static population synthesis 
technique was more popular. 

However, there has been remarkable progress in the advanced stellar evolution theory during 



the past few decades, especially for the core helium-burning phase (e.g. Rood 1973 ; Sweigart, 
Mengel, & Demarque 1974; [Sweigart 1987 ; Lee &: Demarque 199C| ; Castellani &: Tornambe 1991 ; 



Alternatively, "optimized (or optimizing) synthesis" (3'Connell 1987) 
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Horch et aTTggg ; Dorman, Rood, & O'Connell 1993; |Fagotto et al. 1994| ; YDK). Therefore, many 



groups are now including advanced phases of stellar evolution in EPS (e.g. |Nesci &: Perola 1985 ; 



Magris & Bruzual 1993; Bressan, Chiosi, & Fagotto 1994; Buzzoni 1995| ; Dorman, O'Connell, &; 



Rood 1995; Yi et al. 1995; Park Lee 1997| ). Many of them are mainly aiming at solving the 



"UV upturn problem", where the UV upturn is defined as the increasing flux in the UV (1,000 
2,500 A) with decreasing wavelength, as found in giant elliptical galaxies ( |Code fc Welch 197S ; 
Faber 198$ ; |Burstein et al. 1988Q . 



The most serious concern about such elaborate population studies is that the effects of input 
parameters on the resulting UV flux are very poorly understood. Unlike population synthesis in 
the optical band, UV population synthesis for old stellar systems, such as elliptical galaxies, is very 
sensitive to various detailed physical assumptions. This is because most candidates for UV sources 
are highly evolved stars whose properties are governed by several input parameters that have not 
yet been determined very well. Among such parameters are mass loss and its distribution, galactic 
helium enrichment (AY/AZ), and initial mass function. In particular, as we will show later, the 
assumptions about mass loss and its distribution influence the magnitude of the resulting UV 
flux significantly. Many recent population synthesis studies have been presented to the readers 
without addressing the sensitivity of their results and conclusions to alterations in the adopted 
input parameters. 

The main goals of this paper are to illustrate the sensitivity of the model UV spectrum to 
various input parameters, and to investigate the plausibility of UV-bright core helium-burning 
stars as the major UV sources in giant elliptical galaxies. This paper complements previous 
works in the literature. Greggio & Renzini (1990) explored the sensitivity of the UV modeling to 
AY/AZand to the metallicity-dependence of mass loss. The effects of age and metallicity have 
been investigated by several groups, among which the works of Bressan et al. (1994), Dorman et 
al. (1995), and Tantalo et al. (1996) are noteworthy. Although these groups and others (most 



notably the empirical works done by the HUT team; Ferguson & Davidsen 1993; Brown et al.| 



1997) have tested the plausibility of core helium-burning stars as the main UV sources in giant 
elliptical galaxies before, we found it necessary to investigate it within a fuller parameter space, in 
those cases where several input parameters' influence are still poorly known. 

Realistic synthetic HB, based on gaussian mass dispersion on the HB, has been used in our 
population synthesis. Although Park & Lee (1997) used the same technique in their recent study 
for the first time, ours is the first to combine the mass loss as a function of age and metallicity 
and gaussian mass dispersion. We will show why realistic HB treatment is important. We also 
would like to demonstrate that intricate astrophysics must be considered in estimating mass loss 
and in incorporating the estimates into the population synthesis code, even for the most simplistic 
approach (e.g. Refiners' mass loss formula). 

In order to distinguish the effects of these input parameters from others, all the models in this 
paper are based on simple model populations in which all stars are born at the same time with 
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the same chemical composition. This study focuses on the sensitivity of the UV modeling to the 
adopted input parameters, and, thus, does not include a detailed comparison between models and 
observational data. The following paper (Yi, Demarque, &i Oemler 1997) in the series is mainly 
dedicated to such comparisons. 



2. Stellar Spectral Library Construction 

An obvious source of uncertainty in EPS is the stellar spectral library used to convert the 
theoretical quantities (such as T e g , log g , and radius R) into observable ones (e.g. magnitudes, 
colors and spectral line strengths) (Chariot, Worthey, & Bressan 1996). Galaxy models have 
been based on either empirical stellar spectral libraries which are highly incomplete in the super 
metal-rich (Z > Zq ) regime and lack UV data, or on theoretical libraries that do not reproduce 
observed stellar spectra precisely, and have not been fully verified empirically at UV wavelengths. 



There are several comprehensive empirical libraries (e.g. |Gunn fc Stryker 1983 ; Pickles 1985a ; 



Silva fc Cornell 1992| ) in the optical through near-infrared (IR) bands, and they have been used 



in many population synthesis studies. In general, they are adequate for MS through RGB stars 
and in the solar abundance regime. However, elliptical galaxies are believed to be metal-rich, and 
no adequate empirical library exists for super metal-rich {Z > Zq ) stars. In particular, the lack 
of empirical spectra of evolved, metal-rich stars, such as SBP (slow blue phase) stars, is a serious 
drawback for UV studies. The best UV spectral library is derived from IUE observations (Fanelli, 
O'Connell, & Thuan 1987; |Fanelli et al. 1992| ); however, it lacks super metal-rich (> Zq ) samples, 



and does not cover the far UV near the Lyman limit that is crucial to the UV upturn study. 

On the other hand, theoretical spectral libraries suffer from insufficient atomic data and 
physical understanding. Even the most recent models may be inadequate, especially for cool stars 
( Buser &: Kurucz 1992 ; Morossi et al. 1993 ; Kirkpatrick et al. 1993| ; Gustafsson fc J0rgensen 



1994 ). These theoretical libraries have to be validated by acquiring as much empirical data as 
possible. But, the validation in the UV wavelength is possible only using the space facilities which 
have limited and competitive observing time. 

While both empirical and theoretical libraries have problems, a great advantage of using the 
theoretical spectral library is that we can clearly see the effects of various parameters, such as 
temperature, metallicity, and surface gravity, in the models. Moreover, the recent Kurucz library 
covers the high metallicities and high temperatures that are crucial to highly evolved stars and 
thus to UV studies, but that are not available empirically. Thus, a theoretical library has been 
chosen, assuming that the uncertainty in the model spectra is still smaller than the effects of 
temperature and metallicity that cannot be studied adequately using empirical libraries. 

The foundation of the spectral library in this study is the Kurucz library (1992). It covers 
[Fe/H] = -4 - +1, T cS = 3,500 - 50,000 K, and log g= 0.0 - 5.0 (see Kurucz (1992) for details 
of the grids). Additional spectra have been constructed for the stars with high T e g and/or high 
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log <? using the 1995 version of the Hubeny spectral synthesis code, TLUSTY178 (see Hubeny 1988 
for details). The effects of T e g- and metallicity are demonstrated in Figures 1-2. 

Since the metallicity effects are small when T e g ^ 50,000 K, all models for T c q > 50,000 K 
have been constructed for Z = Zq . The stellar spectrum for each star is constructed using linear 
interpolation in metallicity and surface gravity. However, for T e g , the population synthesis code 
finds the model spectrum with the closest temperature in the library instead of interpolating in 
the temperature grid in order to save CPU time. For example, the code will use model spectra 
of T e fj = 3,500 K and 3,750 K for the stars of T e g- = 3,400 K and 3,650 K. This approach slightly 
overestimates the contribution from such cool stars in the composite spectrum in the IR because 
the minimum T e g- in the Kurucz library is 3,500 K. However, this effect is negligible in the UV and 
optical range because only a small number of stars are that cool and such cool stars contribute 
little to the integrated light in the UV to optical range. 



3. Construction of synthetic CMDs and SEDs 

EPS consists of two steps, constructing model color-magnitude diagrams (CMDs) and spectral 
energy distributions (SEDs). It will be shown later in this study that a careful synthetic CMD 
construction including a realistic HB is crucial, especially for UV population synthesis. The model 
CMD provides basic information, such as [Fe/H] , T c g , \ogg , and the radius of the star R. Since 
a model stellar spectrum is defined only by [Fe/H] , T e g , log g , and R, i.e. 

fx = f([Fe/H], T eff , log g, R), (1) 

the convolution of this information with the spectral library generates a model composite 
spectrum. 

The core helium-burning phase in the population synthesis requires a complicated treatment 
involving a synthetic HB construction. Therefore the core hydrogen-burning phase (MS and RGB) 
and the core helium-burning phase - strictly speaking, post-RGB: HB, asymptotic giant branch 
(AGB), slow blue phase (SBP), post-AGB (PAGB), and white dwarf (WD) - have been dealt with 
separately. Then, the total flux F\ is the sum of the flux from the core hydrogen-burning stars F^ 
and that from the core helium-burning stars F^ e 

F X = Jf + Ff e . (2) 

The foundations of the model CMD construction are the Yale Isochrones 1996 ( pemarque et 
al. 1996| ) and Yi et al.'s post-RGB tracks (YDK). 
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3.1. Core hydrogen-burning phase: MS and RGB 

The core hydrogen-burning phase is relatively easier to deal with because both MS and RGB 
are well understood. The Yale Isochrones 1996 provides all the necessary information of the 
core-hydrogen burning stars for the population synthesis. 

The initial mass function (IMF) is assumed to follow a conventional power law as follows, 

^^M"<-) (3) 

where Mi is the stellar mass, rii is the number of stars in the ith mass bin, a is a normalization 
constant, and x is the IMF slope. Three values of x, namely -1, 1.35, and 3, are chosen to see the 
effects of IMF on the model CMDs and SEDs; the models with x = 1.35 ( |Salpeter 1955 ; Larson 
19921 ; |Larson 19951 ) are considered to be standard. 



The flux from all the core hydrogen-burning stars, on the MS and RGB, , is 

Fx =5>i fx,i (4) 

i 

where 

i-Mi+SM 

n i = a M-^ +1 ) dM. (5) 

JMi-5M 

The normalization constant a has been set to make the mass of the model galaxy 10 12 Mq . 
This large number avoids any stochastic effect from rare bright stars. One of the nice features of 
the population synthesis of giant elliptical galaxies is that giant elliptical galaxies are composed of 
a sufficiently large number of stars that we can avoid a serious stochastic effect where statistical 
fluctuations of a few bright stars can influence the SED. On the contrary, population synthesis for 
small galaxies and star clusters suffers from the stochastic effect. Therefore, such studies must be 
approached with caution. 



3.2. Mass loss estimation 

Mass loss is one of the most important, but least understood processes in stellar astrophysics. 
It has a dominant influence on the UV evolution of galaxies because the mass loss determines the 
mean T e g of core helium-burning stars that are the major UV sources in old stellar systems, as 
will be shown later. 



Reimers' empirical formula of mass loss ( Reimers 1975| ) 



dM „ L 

_ = _4xl0-%_, (0) 

where L is the luminosity, g is the surface gravity, and R is the radius, has been used in this study. 
The mass loss efficiency parameter, 77, an empirical fitting factor, was originally chosen to be unity 
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based on the study of 16 metal-rich (Population I) red giant stars (Reimers 1975). However, r\ has 



been reported to vary remarkably from 0.25 to 2 - 3 ( pupree 1986| ; Kudritzki & Reimers 1978 



Renzini 1981). While the estimated r\ for metal-rich stars varies significantly, the estimated range 
of r] for metal-poor stars seems to be much better determined. Most studies on metal-poor stars 
suggest rj = 0.3 - 0.6 ( [Aaronson & Mould 1982] [Mould fc Aaronson 1982] ; [Renzini 198~l[ ; Lee, 
Demarque, & Zinn 1994). 

One of the best ways to estimate the mass loss in metal-poor stars uses the fact that the mass 
of the RR Lyrae variables (a type of HB star) and the mass of the red giants can be determined 
independently. In principle, the mass loss on the RGB can be deduced from the mass difference 
between red giants and HB stars. However, a precise, empirical mass estimation of ordinary HB 
stars is difficult to achieve. Several pulsation studies showed that the mass of a particular type of 
RR Lyrae stars (RR-d double mode variables) can be reasonably estimated ( Petersen 1973; ) . The 
mass of RR Lyrae stars can also be measured via stellar evolution theory by fitting the observed 
HB morphology with synthetic HBs (Yi, Lee, & Demarque 1993). For a long time, there existed 
an inconsistency between the evolutionary mass of RR Lyrae variables and their pulsation mass, 
in the sense that the pulsation masses were lower (Petersen 1973; Lee, Demarque & Zinn 1990). 
This discrepancy raised doubts in the minds of many variable star researchers about the validity of 
the HB evolutionary models (e.g. Clement, Kinman, & Suntzeff 1991). Accordingly, the true mass 
loss on the RGB was unclear. However, the mass estimates from the two independent methods 
are now in good agreement when improved physics including the OPAL opacities is used (Kovacs, 
Buchler, & Marom 1991; |Cox 1991| ; |Yi et al. 199^ )Fl Because of this agreement, we can reasonably 



estimate the amount of mass loss experienced by the red giant progenitors of double mode RR 
Lyrae stars, once the masses of red giants are provided. The masses of the red giants as a function 
of age and chemical composition can be reliably deduced from the stellar evolution theory. 

In order to redetermine r\ for metal-poor stars in this study, we adopted the evolutionary mass 
of RR Lyrae stars in metal-poor globular clusters, M15 (Z = 0.0001) and M3 (Z = 0.0004), from 
Yi et al. (1993) and the mass of red giants from the Yale Isochrones 1996. It should be emphasized 
that rj estimation based on the HB morphology, that is, on the mass estimates of RR Lyrae stars, 
must use the new mass estimates of RR Lyrae stars after the RR Lyrae mass discrepancy has been 
removed. If the old, lower mass is used, the r\ would be overestimated by almost a factor of two. 

The mass loss (Table 1-4) has been calculated using Reimers' formula which is applied 
to MS - RGB stellar evolutionary tracks that have been constructed without mass loss. Several 
evolutionary tracks of different masses were selected for the mass loss estimation in order to 



4 There still may be inconsistency in the estimated masses of HB stars. Recent spectral analyses (e.g. Moehler, 
Heber, & de Boer 1995; de Boer, Tucholke, & Schmidt 1997) of hot stars on the blue tail of M 15 report that the 
mass estimation is substantially lower than what canonical HB models suggest, raising a question either about the 
origin of the hot stars on the blue tail or the validity of stellar evolutionary calculations. On the other hand, de Boer 
et al. note that "the likely cause for the discrepancy of the mass values lies in the low gravities derived from Balmer 



line profiles." (see also Cacciari et al. 1995) 



- 8 - 



approximately cover the range of 1 - 25 Gyr of age. For example, in the case of Z = 0.02 (~ Zq ), 
the tracks of 1.8 M Q (lifetime on the MS - RGB « 1.3 Gyr), 1.4 M Q (3.3 Gyr), 1.0 M Q (11.1 Gyr), 
and 0.8 M Q (31.5 Gyr) were used to estimate the mass loss at 1.3, 3.3, 11.1, and 31.5 Gyr, and 
linear interpolation was carried out within this age range. Therefore, in a 10 billion years old 
stellar system with (Z, Y) = (0.02, 0.27) where most stars on the MS through RGB are slightly 
heavier than 1.0M & , Reimers' formula suggests that stars experience mass loss of about 0.119 
(if n = 0.3) - 0.396 Mq (if r/ = 1.0) (see the second top left panel in Figure 4 and Table 1-4) 
depending on the adopted rj. We assume the mass loss becomes noticeable when L > Lq . Also, 
the mass loss has been forced to stop in the numerical calculations either when the star reaches the 
helium core flash stage or when it causes the envelope mass M env to be smaller than 0.01 Mq (as 
the star loses mass and as the core mass M core grows). 

It has been suggested that M env cannot be infinitesimally small ( [Swcigart et al. 1974| ) for 
the following reason. A certain mass is required to be present in the stellar envelope in order 
to exert sufficient gravitational pressure on the degenerate core of an RGB star and to cause a 
helium core flash (we call this the minimum M env hypothesis). This is how HB stars are born. It 
is believed that nearly all the energy produced by the helium core flash is consumed in removing 
the degeneracy of the core and that there is not enough energy to expel any significant amount of 
envelope material into space ( [Tomasko 197C ; Cole &: Deupree 1980| ; Cole, Demarque, & Deupree 



1985). The minimum M env has been suggested to be approximately 0.01 - 0.02 M & (Hayashi, 
Hoshi, & Sugimoto 1962; |Refsdal fc Weigert 1970| ; [Sweigart et al. 1974j) . 



However, it is still debatable whether such a minimum M env exists considering that there are 



hot (low-mass) field subdwarf stars (see Sweigart et al. 1974 for review) and low-mass UV-bright 



stars in NGC6791 flKaluzny fc Udalski 199^ ; Liebert, Saffer, & Green 1994; [Kaluzny & Rucinskj 



1995) that seem to have an envelope of M env <^ 0.01 Mq . In particular, D'Cruz et al. (1996) 
recently suggested that such a minimum M env does not exist. Instead, they claimed that red 
giants with a very smaller envelope mass experience helium-core flash after leaving the tip of the 
RGB quietly. We do not include such models in this paper, but we believe that their finding 
is quite plausible and deserves further investigations both to verify it with thorough theoretical 
calculations and to explore the effect of it to the population synthesis. While the value of the 
minimum M env is not yet clear, 0.01 Mq has been chosen in this study in order to approximately 
accommodate the hot stars in NGC 6791 and the hot field subdwarf stars. 

Much of the mass loss is believed to occur during the bright phase of the RGB because it is 
sensitive to the stellar radius R, and R increases dramatically as a red giant evolves along the 
RGB. If the minimum M env hypothesis is removed, a larger mass loss would be allowed for some 
stars and the synthetic HB would contain hotter stars. However, the absence of such hot HB stars 
in Galactic globular clusters in general supports the minimum M env hypothesis^]. Table 5 lists the 



5 While there are several clusters that are known to contain quite hot HB stars (e.g. M 15, NGC 6752, NGC 1904, 
& uj Cen), only a few of such hot stars are believed to have an M env smaller than 0.01 Mq . If such hot stars with 
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estimated mass of the RR Lyrae variables Mrr, the mass of the RGB stars at the tip Mrg, an d 
the best fitting r\ for the metal-poor stars. 

The estimated rj depends on the assumed ages of the globular clusters used in this study. We 
have chosen 12 - 16 Gyr as an acceptable age range of the Galactic globular clusters M15 and M3 
following Chaboyer, Demarque, & Sarajedini (1996). The line with asterisks in the top panels of 
Figure 3 is the locus of the mass loss in the adopted age range, and this yields the estimated mass 
loss efficiency of r/ ~ 0.2 - 0.5, which is in good agreement with the estimates from the previous 
studies of metal-poor stars as introduced above. Previous studies (see above) seem to suggest 
that estimated ry's for metal-poor stars are more consistent with one another than for metal-rich 
stars, in which r/ varies from 0.25 to 2 - 3. Despite the large uncertainty, this seems to suggest a 
metallicity-dependence of r\ in the sense of increasing r/ with increasing metallicity if, as suggested 
by Reimers (1975), r] m 1 for Z = Z & . 

Figure 4 shows the mass loss estimates, AM , for metal-rich stars. There are two general 
trends to be noted in Figures 3-4. Firstly, AM generally increases as age increases mainly because 
Mrg decreases with increasing age and the lifetime of a red giant (t#) becomes significantly longer 
as its mass {Mrg) decreases. Secondly, AM reaches a maximum when Mrq becomes fairly low. 
This is because, below this mass, the stellar evolutionary pace is slower than the stellar core 
growing pace. The mass loss is a product of the two competing time scales, the evolutionary 
time scale on the RGB and the core growth time scale. Figure 5 shows that AM increases as 
Mug decreases until a certain mass (Mrg ~ 1 Mq in case of Z = Z & ) is reached because tn 
increases as Mrg decreases and AM obviously increases as tn increases. Below this mass, however, 
the evolutionary pace is much slower than the core growing pace, so before the star reaches the tip 
of the RGB, its core becomes large enough to cause M env < 0.01 Mq . Thus a further evolution on 
the RGB is halted and the post-RGB phase begins. The post-RGB stars with Menv < 0.01 Mq are 
assumed not to become HB stars but to become WD stars following the M = 0.453 Mq track 
of Sweigart et al. (1974). As discussed above, it is possible that some red giant stars with total 
masses very nearly equal to the helium ignition core mass will start core helium burning while or 



after crossing from red to blue giant in the HR-diagram ( Demarque fc Mengel 1971 ; Castellani fc 



Castellani 1995 ; D'Cruz et al. 1996 ). In this crucial mass range, the mass loss estimation is highly 
sensitive to how the mass loss takes place (which is poorly known), and thus it has to be carried 
out with much caution. The same phenomenon exists even if the minimum M env hypothesis is 
removed. 

One can avoid such complications by including mass loss in the stellar evolution calculations 
from the beginning, if the mass loss mechanism were better understood. However, the current 



Menv < 0.01 Mq are remnants of the single stellar evolution, they may provide counterexamples to the minimum 
Menv hypothesis. However, very low mass envelope HB stars might instead be created by mass transfer in binary 
systems (Mengel, Norris & Gross 1976), or by close interactions in dense parts of globular clusters ( |Bailyn 1995 



Fusi 



Pecci et al. 1993). In this paper we have considered only the evolution of single, non-interacting stars. 
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standard models do not yet take into account mass loss simply because the understanding about 
mass loss is still primitive. Reimers' formula suggests most mass loss takes place near the end 
of the RGB where surface gravity becomes smallest. Moreover, theoretical hydrodynamical 
models (Bowen Willson 1991; Willson, Bowen, &: Struck 1996) suggest that mass loss on the 



RGB is an even more abrupt phenomenon near the end of the RGB than predicted by Reimers' 
formula. Willson (1996, private communication) even suggested that it would be closer to truth 
to remove a certain mass from the star at the tip of the RGB than to include Reimers' formula 
in the stellar evolutionary calculations. If this is true, whether a mass loss formula is included in 
the construction of the tracks or not causes little difference on the stellar evolution and on the 
estimated mass loss. It may even be prudent not to include an arbitrary mass loss formula in the 
evolutionary calculations untill a more reliable mass loss formula becomes available. 

The theoretical hydrodynamical models suggest that mass loss takes place when the star 
reaches a certain critical luminosity (the "cliff"), which is a function of mass and metallicity. 
Below the cliff, the mass loss rate is too low to be observable; and above the cliff mass loss takes 
place so rapidly that the star is unlikely to be observed in this phase of evolution. It is significant 
that their models precisely reproduce Reimers's formula with r\ = 1 for Z = Z & . Moreover, they 
explain the positive metallicity-dependence of rj which observations may be indicating. While their 
models have not yet been parameterized for population synthesis, Willson et al. (1996, private 
communication) advise using Reimers' formula with variable r\ in the sense that r\ increases as 
metallicity increases. We call this the variable r/ hypothesis. Table 6 shows the suggested r/ for 
different metallicity. 

The population synthesis models in this study have all been constructed using values of rj = 
0.3, 0.5, 0.7, and 1.0 that roughly span the range of empirical and hydrodynamical estimates. This 
will allow us to investigate the effects of mass loss efficiency on the model CMDs and SEDs. The 
impact of the variable rj hypothesis, which has been investigated earlier (Greggio fc Renzini 1990| ) 



for single abundance models, will be revisited for the composite galaxy models in our next paper 
(Yi, Demarque, & Oemler 1997). 



3.3. Core helium-burning phase: post-RGB 

We adopt the evolutionary tracks of core helium-burning phase mainly from YDK. Since 
YDK do not provide post-AGB (PAGB) tracks, the PAGB tracks of Kiel group ( [Schonberner 



1979| ; tSchonb erner 1983] ; [Blocker &: Schonberner 1990|) have been adopted. The YDK HB tracks 



in this study are the short version (plus symbols in Figures 1 - 2 of YDK) of the full tracks (as 
solid lines in their Figures 1-2). The time step on each HB track is 5 Myr and the mass step is 
approximately 0.04 M & , with smaller steps near the blue end of the HB as listed in Table 1 of 
YDK. The masses of the PAGB stars and WD's and the process that determines them are yet to 
be understood. Meanwhile, Kiel group's models of M = 0.546, 0.565, and 0.598 Mq have been 
used in the actual population synthesis. One of these three models was chosen as the progeny 
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of an HB star (almost arbitrarily) depending on the mass of the HB star; that is, 0.546, 0.565, 
and 0.598 M models were chosen for the HB stars with M < 0.56, M = 0.60, and M > 0.64. 
However, the choice of the model does not affect the resulting model integrated spectrum much 
because PAGB stars are not likely to be the major UV sources when a galaxy shows a strong UV 
upturn, as will be shown in Section 4.1. 

YDK (Figures 1 - 2 of YDK) show that late stellar evolution in the core helium-burning 
phase, especially in metal-rich galaxies, is very complicated. And, such evolved, metal-rich 
core helium-burning stars are potentially significant UV sources. Therefore, in order to fully 
incorporate the complete evolution of core helium-burning stars in the population synthesis code, 
a physically plausible treatment of the phase is crucial. 



The modified gaussian mass distribution ( Lee et al. 1990 ), first invented by Rood (1973) 
and later modified and elaborated on by Lee et al., has been assumed in this study for the core 
helium-burning phase as follows: 



$(M) = V Q [M - (M HB - AM)](M RG - M) exp[- 



(M HB - MY 



(7) 



where a (= 2a s d) is a mass dispersion factor, is a normalization factor, and Mhb (= Mro — 
AM) is the mean mass of HB stars. Three values of a (0.01, 0.06, and 0.99 in M ) have been 
chosen to represent distributions similar to a delta-function, a globular-cluster-function ( [Lee et aT 



1990), and a near-constant-function, respectively. While the true value of a is still quite uncertain 



([Rood 1973| ; |Smith k Norris 1983| ; [Lee et al. 1990] [Walker 1992] |Catelan 1993| ), we have chosen a 
= 0.06 Mq as the "standard" value in this study following empirical (but, with no physical basis) 
choices of Rood (1973) and Lee et al. (1990). The shape of the functions are shown in Figure 6. 
In actual calculations, we first estimate the mass loss as a function of age and metallicity, then, 
apply the gaussian distribution on the mass loss^J. Assumption of a (modified) gaussian function 
for the HB mass distribution has been questioned ( Smith fc Norris 1983| ; [Rood 19*90 ; Ferraro et 
al. 1992| ; pixon et al. 1996|) . However, the single-gaussian distribution has been very successful in 



reproducing the observed HB morphology in Galactic globular clusters ( Lee et al. 1990) ; Lee et al. 
1991 |Dixon et aI.~l99T| ; |Catelan 1993[ ). 



The total number of core helium-burning stars, Nn e , is as follows: 

Nne = tHe RHe 



(8) 



where tHe is the lifetime of the core helium-burning star of M = Mhb , and Rhc is the rate at 
which stars leave the RGB for the core helium-burning stage. So, at a fixed time, there will be 



6 Recently, J0rgensen & Thejll (1993a) and D'Cruz et al. (1996) attempted to explain the width of the HB of 
Galactic globular clusters by means of a spread in r\ (mass loss efficiency) rather than a spread in mass loss. J0rgensen 
& Thejll (1993a) noted that "star-to-star variations of r], between 77 = 0.0 and r/ = 0.7 can explain the HB morphology 
in typical globular clusters" . It would be very interesting to investigate the effect of such HB treatments to the UV 
population synthesis. 
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NH e core helium-burning stars on all post-RGB phases. The rate, Rm in number of stars per unit 
time interval t\ — ti is defined to be 

Rue = / M-^ 1 ) dM. (9) 

h — l 2 JM GB {t 2 ) 

The calculated probability, VE'(M), is applied to each evolutionary track according to its mass 
Mi. There is a different list of stars for different chemical compositions, and each track has a 
number of evolutionary points with evolutionary time scale information^. Thus, the number of 
stars at the jth. evolutionary point on the ith track will be 

mj = *(M) N He ^ (10) 

and 

U = E *iJ ( n ) 

3 

where Vl/(Mj) is the probability of the star of M = Mi to be present, ti j is the time spent by the 
ith star in the table near the jth position in the CMD and ti is the lifetime of the ith star in the 
post-RGB phase. The sum of the probability ^(Mj) should be unity: 

* = ]T * (Mi) = 1. (12) 



Then the total flux from the core helium-burning stars, F^ e , is 

Fx e = EE fw "iJ. ( 13 ) 

i j 

where fx,ij is the H ux from the ith star in the list of post-RGB stars near the jth evolutionary 
point. 



4. Effects of input parameters 

Spectral evolution of galaxies in the UV is much more uncertain than that in the optical band 
because it is more sensitive to the evolution of post-RGB stars which are less well understood than 
the MS - RGB stars that dominate the optical flux. It is important to understand first the effects 
of input parameters on the resulting CMD and UV SED in order to answer the questions related 
to the spectral evolution of the galaxy, including the UV upturn phenomenon. 

Primary parameters affecting the UV flux are age and metallicity. These are the most 
fundamental quantities in describing stellar populations, because they predominantly determine 



Evolutionary time scale: how long a star stays within a box of given size in the CMD. 
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the T e g 's of post-RGB stars. We would like to investigate whether metal-rich old stellar systems, 
such as giant elliptical galaxies, are capable of exhibiting a pronounced UV flux. Secondary 
parameters include mass loss efficiency, galactic helium enrichment parameter, form of the 
dispersion in mass loss, and the slope of IMF, while these parameters are not necessarily less 
important in producing UV flux than the primary parameters. We assume in this study that all 
stars in a model stellar system form in an instantaneous starburst. 

4.1. Age and metallicity: primary parameters 

The effects of age and metallicity in old stellar systems have recently been widely recognized 
through studies of Greggio & Renzini (1990), Bressan et al. (1994), and of Dorman et al. (1995). 
In this section, we would like to provide a simple explanation of the basic conclusions of such 
previous studies by means of stellar evolution theory. We also illustrate that core helium-burning 
stars are a more likely major UV source in old stellar systems than PAGB stars. Both a higher age 
and a higher metallicity cause a galaxy to look redder in the optical band mainly because MS and 
RGB stars become cooler as they age and as metallicity increases, as is clearly shown in Figures 7 
- 12. The effects of age and metallicity on the UV flux are more complicated, particularly because 
the major UV sources vary with time. Before we discuss the details, it is important to remember 
that we are investigating the effects of age and metallicity while other parameters, such as rj, 
are kept constant. If other parameters are age-dependent and/or metallicity-dependent the true 
effects of age and metallicity become much more complicated. 

4.1.1. Age 

The major UV sources vary with time. As seen in Figures 7 - 12, PAGB stars are the 
only far-UV sources in the early (^5-10 Gyr, depending on the metallicity) history of galaxy 
evolution, when the majority of core helium-burning stars are still too cool to produce any 
UV light. However, once the galaxy is old enough to develop a substantial number of hot core 
helium-burning stars, those stars dominate the UV spectrum. 

In order to investigate the light contribution from various sources, we first arbitrarily divide 
the core helium-burning stage into two: central helium-burning stage and shell helium-burning 
stage^]. The central helium-burning stage is defined as the stage where most of the helium-burning 
takes place near the center of the star. This stage lasts from the zero age HB (ZAHB) to the point 
where the central helium is almost exhausted (Y c = 0.01). Similarly, the shell helium-burning 
stage is defined as the stage that stretches from this point until the star reaches either the tip of 



The terms, "core helium-burning" and "central helium- burning" , may be confusing. Readers should pay attention 
to the different definitions of the terms. 
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the AGB or the beginning of the WD phase, whichever comes first. So the central helium-burning 
stage is generally referred to as the HB, and the shell helium-burning stage as the post-HB phase 
including the evolved HB, the AGB and the slow blue phase (SBP). These two stages are denoted 
as crosses and squares in Figure 13. 

Assuming that all the stars in elliptical galaxies have formed in an instantaneous burst, 
Figure 14 shows that central and shell helium-burning stars are far more efficient UV sources than 
PAGB stars in old galaxies. The upper and lower panels of Figure 14 show the light contribution 
from various evolutionary stages in young and old metal-rich model galaxies. It shows that MS 
and RGB stars are the dominant sources in the visible range at all times. The far-UV spectrum 
is, however, mainly dominated by highly evolved stars. At small ages, PAGB and MS Turn-Off 
stars are the dominant far-UV and near-UV sources (the upper panel), respectively. However, 
such galaxies do not exhibit any significant amount of UV flux, as shown in the 5 Gyr old models 
in Figures 10 - 12. The critical moment which divides the periods in which PAGB stars or core 
helium-burning stars dominate the far-UV flux is a sensitive function of chemical composition and 
other parameters (such as 77). In the case of the extreme- UV (A <J 300A), PAGB stars always 
dominate because of their high T e g 's. However, once galaxies become old enough to contain hot 
HB stars, shell helium-burning stars dominate the flux in the far-UV (A ~ 300 - 1,500 A), while 
both central and shell helium-burning stars are equally important in the near-UV (A ~ 1,500 - 
3,000 A). 

The effect of age on the UV flux is relatively simple. Figures 15 - 16 show the UV-to-V colors 
as a function of age and metallicity. The model magnitudes are defined as m\ = — 2.5 log <f\> 
where < f\> is the mean flux in the bandpass. The < /(1500) >, < /(2500) >, and < /(V) > 
are defined by averaging the flux within the ranges 1250 - 1850 A, 2200 - 2800 AflD orman et al 



1995), and 5055 - 5945 A( |Allen 1976| ), respectively. The observational data are from Table 1 



2 of Dorman et al. (1995). From their table, we selected the clusters for which both m(1500)-V 
and m(2500)-V are available. They are mostly ANS or OAO-2 data. The UIT data of NGC 1904 
(M79) and NGC 5139 (uCen) have also been adapted. As a result, three clusters (NGC 1904, 
NGC 5139, and NGC 5272) have two data points in our adapted cluster sample. We excluded the 
IUE data because the IUE aperture is so small that IUE data do not represent the integrated light 
properly ( Dorman et al. 1995| ). In case of metal-poor systems (Figure 15), UV-to-V colors become 



redder with age for a while as the MS Turn-Off stars, the major UV sources in young metal-poor 
systems, become redder with time, and then they turn around to become bluer as HB stars get 
hotter. Most metal-poor globular clusters are reasonably matched by the models with rj = 0.5 
but metal-rich globular clusters like NGC 6388 and 47 Tuc are not matched by the same models^]. 
This issue of r/ will be discussed in the next section. Similarly, metal-rich elliptical galaxies are 
reasonably reproduced by the models with ry = 0.7 in Figure 16. 



9 The m(1500)-V of 47 Tuc obtained by OAO-2 is likely to be the lower (bluest) limit because there is one bright 
PAGB star in OAO-2's large aperture (Dorman et al. 1995). 
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4-1.2. Metallicity 

The contribution from each evolutionary stage to the UV flux also varies with metallicity 
because the post-HB evolutionary pattern differs according to metallicity, as shown in Figures 1 - 
2 of YDK. Since the UV-bright SBP phenomenon becomes far more prominent in more metal-rich 
stars, the contribution from the shell helium-burning stars in the far-UV is more conspicuous in 
metal-rich models (Figure 17). 

It should be noted that the light contribution from evolved core helium-burning stars to the 
total U band flux (A ~ 3,000 - 4,000 A) is not negligible. The lower panels of Figures 14 and 17 
suggest that substantial amount of the total light comes from such evolved stars. However, note 
that the precise quantity highly depends on the adopted input parameters. Earlier population 
synthesis studies (e.g. [Tinsley 1972j ; p'Connell 1976| ) noticed that population synthesis models, 
that did not include such evolved stars in an adequate way, were predicting much redder U — B 
colors than observed. They suggested that this mismatch could originate from the presence of 
young stars in giant elliptical galaxies. However, Gunn, Stryker, & Tinsley (1981) recognized the 
lack of empirical support for the young star interpretation, and suggested blue stragglers and what 
we now call PAGB stars as possible candidates for the observed UV flux. We understand now, 
as clearly shown in Figures 14 and 17, that a proper treatment of evolved stars, especially with 
a realistic mass dispersion on the HB, enhances the UV - U band flux significantly, as pointed 
out by Burstein et al. (1988) earlier^. This short wavelength region is extremely sensitive not 
only to age and metallicity of the galaxy but also to various input parameters, including mass loss 
efficiency and the HB mass distribution, as will be shown in the following sections. 

The overall effect of metallicity is more complicated than that of age because it is a product of 
several competing effects. Even when we consider only the period in which the core helium-burning 
stars are the main UV sources, it is better to see the metallicity effect in two different metallicity 
domains, (1) Z < Z Q and (2) Z ^ Z Q , 

When Z < Zq , UV-to-U flux ratios decrease with increasing metallicity. This is a result of 
the following competing effects. (1) The temperature of the HB stars decreases dramatically as 



metallicity increases, (2) more metal-rich stars stay on the MS - RGB longer ( Grcggio fc Renzini 



1990 ; [J0rgensen fc Thejll 1993b| ), but (3) a larger mass loss of more metal-rich stars for a fixed rj 



(Figures 3-4) compensates for these two effects ( |j0rgensen fc Thejll 1993b ). The phenomenon 



(2) is mainly caused by the delay of the evolution due to the opacity effect. Consequently, the 
phenomenon (2) makes a more metal-rich cluster have more massive red giants at a fixed age, 
but (3) compensates for this effect and makes the Mhb of more metal-rich clusters smaller than 
that of metal-poor clsuters in this metallicity domain. However, the phenomenon in (1) again 
negates the result of the competition between that of (2) and (3). The net effect of all these 



10 The referee pointed out that the "excess" near-UV (U band) light in elliptical galaxies could also be due to 
intermediate temperature (younger) main sequence stars. 
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competitions is a decrease in the UV-to-V flux ratios with increasing metallicity which is consistent 
with observations, as illustrated in Figure 15. In Figure 15, three metal-poor models covering 
1-25 Gyr of age are compared with globular cluster data. It is clear that both observations 
and models indicate that more metal-rich systems in this metal-poor regime are redder in these 
UV-to-V colors. It should be noted that this trend (the correlation between the UV-to-V ratios 
and metallicity) is based on the assumption that other parameters are identical. In the case of 
Galactic globular clusters, it has been suggested that there is a substantial spread in age among 
them which influences the HB morphology and thus the flux ratios. A large scatter in the flux 
ratios in Figure 16 may have been caused by such differences in cluster parameters. 

In the metal-rich regime Z ^ Zq , the evolutionary phenomena are quite different. Firstly, 
the higher the metallicity, the faster the MS - RGB evolution, because the luminosity from the 
hydrogen burning is significantly sensitive to the mean molecular weight that is determined by the 
chemical composition. This has been noticed previously by a number of workers (e.g. |Greggio~& 



Renzini 1990| ; [Igr gensen fc Thejll 1993b|). Under the assumption of a positive AY / AZ (~ 2 - 3 in 



this study), a higher metallicity means a higher helium abundance. A higher helium abundance 
causes a higher mean molecular weight that results in a faster stellar evolution (see YDK). This 
causes the Mug of more metal-rich stars at a fixed age to be smaller. Secondly, the increase in 
mass loss as a function of metallicity (Figures 3-4) makes Mhb of more metal-rich HB stars much 
smaller at a fixed age. In addition to this, even many of the metal-rich HB stars that are born 
as extremely cool stars quickly evolve to become UV-bright stars (namely, SBP stars in Figures 
1 - 2 of YDK). As a result, more metal-rich galaxies are expected to contain more UV-bright 
stars at a given age in this metal-rich regime, given other parameters fixed. This is illustrated 
in Figure 16: the more metal-rich, the bluer. Figure 16 seems to suggest that the galaxy data 
are better fitted by the Z = 0.01 model than more metal-rich models. We believe that this is 
an artifact that is produced by the comparison between the galaxy data and single-abundance 
models. Obviously, galaxies are not composed of stars of a single abundance only. The effect of 
the composite population treatment will be discussed in detail in our next paper. 

The data and models in Figure 18 show that metal-poor globular clusters in general exhibit 
bluer m(1500)-V and m(2500)-V colors than elliptical galaxies (bottom panel of Figure 18). This 
is mainly because of the opacity effect on the HB stars (in the case of m(1500)-V) and on the 
MS stars (in the case of m(2500)-V). For a given age, the opacity effect on the stellar evolution 
causes metal-poor blue HB stars to have higher T e g- 's than metal-rich counterparts, even though 
the mean mass of metal-rich HB stars, Mhb , are likely to be smaller (note that, given the other 
parameters fixed, a lower-mass HB star is hotter) according to Refiners' formula (see Figures 
3-4 and Table 1 -4). Moreover, even if their T e g 's are identical, the opacity effect (via the 
line-blanketing effect in the spectral synthesis) causes more metal-rich HB stars to look much 
redder than metal-poor HB stars. For the same reason, metal-poor systems show bluer m(2500)-V 



colors, but this time because of the opacity effect on MS stars ( Dorman et al. 1995 ). As a 



result, metal-poor populations exhibit stronger UV-to-V flux ratios than metal-rich populations. 
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However, metal-poor populations with such blue UV-to-V colors do not exhibit a "UV upturn" 
(the continuously increasing flux with decreasing wavelength below 2,500 A). Instead, their UV 
spectra, even for the UV-strong globular clusters, are rather flat below 2,500 A. Consequently, 
globular clusters have redder m(1500)-m(2500) colors than elliptical galaxies, as shown in the top 
panel of Figure 18. So, it should be said that only elliptical galaxies exhibit a "UV upturn" in a 
strict sense, as models predict in Figure 18. This result disagrees with Park & Lee (1997) 's recent 
suggestion. Details will be discussed in our next paper. 



4.2. Mass loss efficiency 

Figures 19 - 20 show that the magnitude of the UV flux is sensitive to the mass loss efficiency 
n in Refiners' formula. A larger n results in lower-mass, hotter HB stars. Therefore, a value for 
mass loss efficiency should be adopted carefully and with sufficient justification. The probable 
range of n seems to be 0.3 -1.0 (or even 2) with a possible positive metallicity-dependence 
supported both by observations and hydrodynamical models, as discussed in Section 3.2. Such 
sensitivity of the UV flux to the mass loss efficiency has significant implications. 

If Reimers' formula is an adequate approximation of reality, n of stars of Z ;> Z & should 
be near or larger than 0.7 (this value depends on the assumed metallicity of giant elliptical 
galaxies) in order for a metal-rich elliptical galaxy to exhibit a UV upturn within a Hubble 
time^]. Old, metal-rich galaxies should have an appreciable number of UV-bright stars and a 
noticeable UV upturn as long as n ^ 0.7 for metal-rich [Z ^ Z & ) stars. On the other hand, if 
n 0.7 for metal-poor stars, as both observations and theory suggest, metal-poor stars cannot be 
major UV sources in elliptical galaxies because they do not produce sufficient far-UV flux (e.g. 
m(1500)-m(2500) color) within a Hubble time, as Figure 20 shows. This problem becomes worse 
when even shorter wavelength regions 1,500 A) are compared, because the UV spectral energy 
distributions of giant elliptical galaxies have much higher far-UV flux than near-UV flux, and this 
cannot be produced by metal-poor models unless n is unrealistically high. Even if larger 77-models 
(e.g. 1.0) may reproduce such high far-UV flux (bottom panel in Figure 20), the predicted overal 
spectral shape in the UV does not look similar to what is observed. Detailed quantitative study 
will be shortly presented by Yi et al. (1997). 

Conversely, the magnitude of the UV upturn can be used to set a constraint on the mass loss 
efficiency of Galactic globular clusters whose age and metallicity are independently estimated^. 
For instance, if we compare the far-UV spectrum of M 79 (Figure 5 in Dixon et al. 199lf ) with 



11 We assume that a Hubble time is slightly larger than the age of the oldest Galactic globular cluster estimated 
using the same stellar evolutionary calculations as the ones used in this study. This study is based on the stellar 



models that suggest an age of approximately 15 Gyr for the oldest Galactic globular clusters (Chaboyer et al. 1996) 
12 In order for this method to be reliable, far-UV (< 1,500 A) spectra should be available. 
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Figure 20, the model with 77 ~ 0.5 reproduces the empirical spectrum best, assuming 15 Gyr for 
its age and "standard" values of model parameters [x = 1.35, a = 0.06). 

Figure 21 illustrates the ^-dependence of metal-poor models. As discussed earlier in Section 
3.2, low-77 models (77 = 0.3 - 0.5) reasonably match the observations at an age of ~ 15 Gyrs. 
However, the metal-rich clusters (filled squares, e.g. NGC 6388) cannot be matched by the models 
with consistent ages. For example, if rj = 0.3 even for such metal-rich globular clusters, UV color 
fitting would suggest NGC 6388 is almost 25 Gyrs old, which is much higher than the average 
isochrone-age of globular clusters (~ 15 Gyr, |Chaboyer et al. 1996 ). If a larger 77 is adopted 



for such metal-rich globular clusters instead, their colors can be reproduced successfully at a 
reasonable age. For instance, the age estimate for NGC 6388 would be ~ 19 Gyr and 15 Gyr 
if 77 = 0.5 and 0.7, respectively. If other complexities are ignored at the moment, this suggests 
77 ~ 0.7 for Z ~ 0.006 (the metallicity of NGC 6388), which is higher than the best fitting 77 for 
metal-poor clusters. This metallicity-dependence of 77 is consistent with the variable-?? hypothesis 
discussed in Section 3.2. Elliptical galaxies are also better matched by high-77 models (77 ^ 0.7), 
as shown in Figure 22, which is again consistent with the variable-77 hypothesis. Note that in this 
diagram we use m(1500)-m(2500) as the y-axis because m(1500)-m(2500) describes the UV upturn 
strength of galaxies better. Regardless of metallicity, 77 certainly plays a very important role in 
model UV flux in the sense that a higher-77, at least in the rage 77 = 0.3 - 1.0, leads to a higher 
UV flux. 



4.3. AY/AZ 



The galactic helium enrichment parameter, AY/AZ , is another important determinant of the 
UV flux, not only because it affects the pace of evolution of MS and RGB stars (|Greggio fc Renzini 



1990 ; l0rgensen fc Thejll 1993b ), but also because metal-rich stars with a higher AY/AZ become 
UV-bright more readily (Horch et al. 1992; Dorman et al. 1993| ; YDK). This phenomenon is 
more significant when Z ;> Z & , with an implicit assumption of a positive AY/AZ . Therefore, a 
metal-rich galaxy with a higher value of Ay / AZ produces a higher UV flux, as shown in Figures 
23 - 24. Qualitatively, models with AY/AZ = 2 develop a similar magnitude of UV flux to that of 
AY I AZ = 3 at approximately 10 - 20 % greater ages. Despite such a difference, it is difficult to 
estimate the true AY j AZTrom UV magnitudes only because the evolutionary paths for different 
Ay / AZ are very similar as shown in Figure 24, unless their ages are already known. 

If AV/AZis not positive (a very unrealistic assumption in the light of what is known about 
nucleosynthesis in massive stars, see e.g. Maeder 1991| ), in other words, if helium does not increase 
with increasing metallicity, several complicated effects will compete with one another: (1) The SBP 
(slow blue phase) phenomenon will not be significant even in metal-rich (but, not unrealistically 
metal-rich) populations ( Dorman et al. 1993| ; Bressan et al. 1994 ). This reduces the UV flux from 
such hot SBP stars. (2) The evolutionary pace of the more metal-rich stars will perhaps become 
slower because the opacity effect will overwhelm the hydrogen luminosity-mean molecular weight 
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effect (See Eqn. 3 in YDK). This is similar to the metallicity effect on the evolutionary pace of 
subsolar-abundance stars; namely, the more metal-rich, the slower the evolutionary pace on the 
MS and RGB if Z <^ Z & . This causes a cooler HB. (3) But, the mass loss will still increase with 
metallicity regardless of helium abundance. This will cause HB stars to become hotter. The final 
result of this competition is not trivial to guess. Unfortunately, the stellar models of AY /AZ = 
are not yet available for numerical experiments. 



4.4. IMF 

The slope of IMF, x, as defined in Eqn. 3, has a very small effect on the resulting UV flux 
as shown in Figures 25 - 27, even if x differs from one galaxy to another significantly. A model 
with a smaller x sends more stars from the MS both to the RGB and to the HB and results in 
an increase in flux in all wavelength regions. So the relative flux is not much affected. Many 
observations have supported x « 1.35 flSalpeter 1955 ; [Miller fc Scalo 197S ; Larson 1992 ; Larson 



1995 ), and thus we have little justification for adopting a significantly different value. The effect 
of the difference in various forms of IMF for the very low-mass stars on the resulting SEDs is 
negligible because they contribute little to the total light of the host galaxy although they could 
affect the calculated mass-to-luminosity ratio significantly. 



4.5. Mass dispersion on the HB 

A smooth mass dispersion on the HB has not been seriously taken into account as an input 
parameter in earlier population synthesis studies^, although it is very important to the resulting 
SED in the UV. Its effect is not monotonic. When a galaxy is young, and has mostly red HB stars, 
allowing a larger dispersion in mass leads to a larger number of hot HB stars, which results in a 
stronger UV flux (Figure 27). On the other hand, when a galaxy is already old enough to have 
most of its HB stars on the hot side, a larger mass dispersion causes a weaker UV flux (Figure 28). 
These effects on the colors are illustrated in Figure 29. As is the case of AY j AZ , it is difficult to 
choose an optimal a from the UV colors because models with even very different values of a follow 
very similar evolutionary paths (Figure 29) unless the age of the galaxies is known a priori. 



13 Dorman et al. (1995) treated the HB as a product of two different values of mass loss on the RGB (i.e. bimodal 
distribution in the mass loss: two peaks - a red HB clump and a blue HB clump - instead of a smooth distribution) 
to understand how large a fraction of hot HB components is needed to account for the observed UV flux of elliptical 
galaxies. Their blue HB clump is either by 5-function or by constant-function within a certain mass range, both 
of which are arbitrary. Although this is a much improved effort compared to previous works in terms of the HB 
treatment, it is still ad hoc because such bimodality is supported (at least so far) only by a few globular clusters, e.g. 
NGC 2808 (Crocker, Rood, & O'Connell 1988). A single- gaussian mass distribution reproduces the HB distribution 



adequately (Lee et al. 199C; Lee et al. 1994), as discussed in Section 3.2 
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Because of this complex sensitivity of the UV flux to the assumed HB mass distribution, a 
realistic HB treatment is essential to the UV population synthesis. Oversimplified HB treatments, 
such as a clump star assumption^, a single-mass assumption^ and a flat distribution] 1 ^, are 
inadequate for UV studies unless there is a sufficient justification. 



5. Summary 

Recently, several groups have been using the evolutionary population synthesis (EPS) 
technique in the hope of understanding the stellar content and the evolution of giant elliptical 
galaxies. Most of these studies are aimed at solving the so-called UV upturn mystery. Although 
EPS in the UV is a powerful technique for these purposes, its dependence on various input 
parameters has not been widely known. Therefore, both modelers and model users are vulnerable 
to misinterpretations of the results from their population synthesis studies. 

The EPS technique has been described in this study in detail. We show that the model UV 
flux is very sensitive to several input parameters as various earlier studies pointed out. These 
input parameters include not only age and metallicity, but also mass loss efficiency 77, galactic 
helium enrichment parameter AY/AZ , and mass distribution on the HB. The effects of these 
parameters are as follows. 

1. Once MS turn-off stars become too cool to dominate the near-UV integrated flux, the 
relative strength of UV flux to visible flux always increases with age because a larger number 
of hotter core helium-burning stars develop as a galaxy ages. Before this point, UV-to-V 
colors become redder with increasing age as MS stars, the major near-UV sources at low ages, 
become cooler. This turning point is sensitive to metallicity and other parameters. It could be 
approximately 1 Gyr or even smaller if Z ;> 0.06 but could be as large as 10 Gyrs if Z = 0.0004. A 
possibility that giant elliptical galaxies have a substantial fraction of very young (< 1 Gyr) stars 
is ignored in this study because of lack of empirical support ( O'Connell et al. 1992| [Bertola et al 



19931) . 



2. In the metal-rich regime (Z ;> Zq ), metallicity has a positive impact on the UV flux 
once a model galaxy is old enough to have a substantial number of hot core helium-burning stars. 
This is caused mainly by the following phenomena in stellar evolution: under the assumption 
of a moderate AY/AZ (~ 2 - 3), (1) more metal-rich stars evolve faster when Z ;> Z & , which 
results in a quicker hot HB development, (2) more metal-rich red giants experience a larger mass 
loss, and (3) more metal-rich core helium-burning stars become UV-bright stars in the slow blue 



14 It assumes that all the HB stars have the same physical properties, luminosities and temperature, ignoring the 
mass dispersion on the HB and the advanced evolution beyond the ZAHB. 

15 It assumes that all the HB stars have the same mass and follow the same evolution. 

16 This assumption allows a mass dispersion but does not consider any concentration of stars at some mass. 
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phase (SBP) more easily. This result is in agreement with Brown et al. (1997) 's recent spectral 
analysis which suggests that "models with supersolar metal abundance and helium best reproduce 
the flux across the entire HUT wavelength range...". Optical spectra of giant elliptical galaxies 
suggest Z ^ Z Q . In the metal-poor regime, the metallicity effect is reversed, however. This is 
mainly because the opacity effect overwhelms the metallicity-mass loss relation. In addition, in 
this metal-poor regime, an increase in metallicity causes stars on the MS - RGB to evolve more 
slowly and causes the HB to look redder. 

3. Because the UV-bright phase (the SBP) of core helium-burning stars is positively 
correlated with helium abundance, a model galaxy with a higher AY / AZ exhibits a stronger 
UV upturn. This agrees with the result of Greggio & Renzini (1990). If, as we expect, the true 
AY / AZ is within the range 2-3, then the degree of the sensitivity seems moderate. The effect of 
AV/AZand the cause of the effect are similar to those of metallicity described above. A model 
with AY I AZ = 2 would develop a UV flux with a similar magnitude to that of a model with 
AY I AZ = 3 at slightly later time (approximately 10 - 20 % of their ages). 

4. A larger mass loss efficiency, 77, causes a stronger UV upturn because it causes a hotter 
HB. The probable range of 77 is 0.3 - 1 (or higher) with a large uncertainty and probably with a 
positive metallicity-dependence suggested by hydrodynamic models. If we admit that the other 
parameters in this study are relatively better known than 77, fixing the other parameters allows us 
to constrain the true value of r\. In the metal-poor regime (Z <J 0.01), if we assume that Galactic 
globular clusters are about 15 billion years old, the IMF slope x « 1.35, and a « 0.06 M & , the 
models with 77 ss 0.3 - 0.5 fit metal-poor cluster data better. However, a larger value-models (e.g. 
77 « 0.5 - 0.7) fit relatively metal-rich globular cluster data (such as NGC6388) substantially 
better. Similarly, under the same assumptions as the metal-poor case, only high-77 (0.7 - 1.0) 
models fit the elliptical galaxy data. 

5. Mass dispersion on the HB is crucial to the UV flux. A larger dispersion causes a stronger 
(weaker) UV flux when the HB stars with the average mass are red (blue). Because the effect 
of a is large, a population synthesis based on a non-realistic synthetic HB is likely to lead us to 
erroneous conclusions. We chose a = 0.06 M from globular cluster studies as our "standard" 
dispersion. However, this value does not have a physical basis and is poorly known for metal-rich 
systems. 

6. Even a significant change in the IMF slope, x, results in only a slight change in the 
resulting UV flux in comparison to optical flux. A smaller value of x causes a slightly stronger 
UV flux, because a model galaxy with a smaller x sends more stars from the RGB to the HB at 
a fixed time. However, there is little justification for an adoption of a much different IMF slope 
from the standard one, i.e. x = 1.35. 
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6. Conclusion 

Based on this sensitivity study, we conclude as follows. Both the positive correlation between 
metallicity and UV-to-V colors among elliptical galaxies ( Burstein et al. 1988| ) and the negative 



correlation between metallicity and UV-to-V colors among globular clusters have been explained. 
Old, metal-rich populations, such as giant elliptical galaxies, naturally develop a UV upturn within 
a reasonable time scale (<^ Hubble time) without the necessity of the presence of young stars. 
Low-mass, core helium-burning (HB and evolved HB) stars are more likely to be the main UV 
sources in the old stellar systems than any other types of stars, as suggested before (e.g. preggic| 
fc Renzini 1990; Horch et al. 1992; Ferguson fc Davidsen 1993j ; [Bressan et al. 199§ ; |D 

orman etj 



al. 1995] ; Yi et al. 1995| ; [Brown et al. 19971 ). For this reason, both metal-poor and metal-rich 



populations can develop a UV upturn if an arbitrarily large age may be assumed. 

Metal-poor models, such as globular clusters, are bluer in m(1500)-V and m(2500)-V than 
metal-rich models like elliptical galaxies as observations suggest, mainly because of the opacity 
effect. But, their UV spectra are flat and fail to exhibit a strong UV upturn because they lack 
very hot HB and post-HB stars. It takes too long (> Hubble time) for a metal-poor population to 
develop a UV upturn of the observed strength if r\ 0.7 for metal-poor stars as both observations 
and hydro dynamical models suggest. For this reason, no metal-poor stellar system has been 
observed to show a strong UV upturn (i.e. the continuous increase in the UV flux with decreasing 
wavelength below 2,500 A) no matter how old it is. Even if metal-poor models with a larger r) 
(v k5 1-0) an unrealistic assumption) may reproduce the observed UV-to-V colors of giant elliptical 
galaxies, the overall spectral shape in the UV would not match the observations very well, as will 
be quantitatively shown in our next paper ( |Yi et al. 1997b| ). On the other hand, metal-rich model 
populations contain a sufficiently large number of UV-bright SBP stars, in addition to hot HB 
stars, to exhibit a UV upturn of the observed magnitude unless the input parameters assumed in 
this study are significantly wrong. Based on this argument, metal-poor stars cannot be the main 
cause of the UV upturn in giant elliptical galaxies that are obviously composite populations, as 
long as the majority of stars in giant elliptical galaxies are metal-rich. 

Two main effects drive a metal-rich population to develop a UV upturn of the observed 
magnitude earlier than a metal-poor population. Firstly, metal-rich stars experience a larger 
amount of mass loss on the RGB when Reimers' formula and a fixed rj are used. This generates 
lower-mass HB stars at a fixed time. This is mainly because a higher opacity causes a larger stellar 
radius. Reimers' mass loss formula (and our intuition as well) suggests that a smaller surface 
gravity caused by the larger radius results in a larger mass loss. Secondly, even if some HB stars 
are not hot on the ZAHB, more metal-rich stars evolve to UV-bright stars more easily. 

In conclusion, it is important to realize that the presence of the UV upturn in the spectrum 
of giant elliptical galaxies becomes neither extraordinary nor unexpected when proper treatments 
of various evolutionary stages are taken into account. It is a very natural consequence of advanced 
stellar evolution that has not been known in detail until recently. While the presence of the UV 
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upturn has been qualitatively understood, the detailed characteristics of the observed UV upturn 
also have to be reproduced in order to confirm the theoretical explanations. All these aspects will 
be investigated quantitatively in the following paper ( |Yi et al. 1997b| ). 

We thank Richard Larson, Robert Zinn, Wayne Landsman, and Sally Heap for useful 
comments. We are grateful to the anonymous referee for constructive criticisms and clarifications 
on many points. This work was part of the Ph.D. study of S.Y. ( |Yi 1996| ) and was supported in 
part by NASA grants NAGW-3563 (S.Y. and A.O.), NAG5-1486 and NAG5-2469 (P.D.). Part 
of this work was performed while S.Y. held a National Research Council- (NASA Goddard Space 
Flight Center) Research Associateship. 



Table 1. Mass loss estimates based on Reimers' formula (77 = 0.3). 



t(Gyr) AM 

(Af ) 

Z 0.0001 0.0004 0.001 0.004 0.01 0.02 0.04 0.06 0.1 0.01 0.02 0.04 0.06 0.1 

Y 0.23 0.23 0.23 0.23 0.25 0.27 0.31 0.35 0.43 0.26 0.29 0.35 0.41 0.53 
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Table 2. Mass loss estimates based on Reimers' formula (rj = 0.5). 



t(Gyr) AM 

(Af ) 

Z 0.0001 0.0004 0.001 0.004 0.01 0.02 0.04 0.06 0.1 0.01 0.02 0.04 0.06 0.1 

Y 0.23 0.23 0.23 0.23 0.25 0.27 0.31 0.35 0.43 0.26 0.29 0.35 0.41 0.53 
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Table 3. Mass loss estimates based on Reimers' formula (rj = 0.7). 



t(Gyr) AM 

(Af ) 

Z 0.0001 0.0004 0.001 0.004 0.01 0.02 0.04 0.06 0.1 0.01 0.02 0.04 0.06 0.1 

Y 0.23 0.23 0.23 0.23 0.25 0.27 0.31 0.35 0.43 0.26 0.29 0.35 0.41 0.53 
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Table 4. Mass loss estimates based on Reimers' formula (rj = 1.0). 



t(Gyr) AM 

(Af ) 

Z 0.0001 0.0004 0.001 0.004 0.01 0.02 0.04 0.06 0.1 0.01 0.02 0.04 0.06 0.1 

Y 0.23 0.23 0.23 0.23 0.25 0.27 0.31 0.35 0.43 0.26 0.29 0.35 0.41 0.53 
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Table 5: Best fitting 77 estimates for metal-poor stars. 



z 


M RR a 


Age(Gyr) b 


M RG c 


AM d 


V 


0.0001 


0.760 


12 


0.852 


0.092 


0.3 - 0.4 


0.0001 


0.760 


14 


0.816 


0.056 


< 0.3 


0.0001 


0.760 


16 


0.788 


0.028 


< 0.3 


0.0004 


0.710 


12 


0.862 


0.152 


0.5 


0.0004 


0.710 


14 


0.826 


0.116 


0.3 - 0.5 


0.0004 


0.710 


16 


0.796 


0.086 


< 0.3 



"Evolutionary mass of RR Lyrae type-d variables from Yi et al. (1993). 

^Because of uncertainties, three ages have been chosen for the clusters, M15 and M3. 

c Mass of red giants at the tip of the RGB, adopted from the Yale Isochrones 1996. 

d AM=M RG ~M RR 
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Table 6: The 77 suggested by hydro-dynamical models. 



Z suggested r\ 

0.0002 0.17 - 0.25 

0.002 0.22 - 0.33 

0.006 0.33 - 0.50 

0.02 1.0 
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Fig. 1. — Variation of SED as a function of temperature for Z = Zq and log 5 = 5.0. From the top, 
T eff = 100,000, 50,000, 30,000, 10,000, 5,000, & 3,500 K. Model spectra of T cS > 50,000 have been 
constructed using Hubeny's spectral synthesis code while others are all from the Kurucz spectral 
library. 
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Fig. 2. — Variation of SED as a function of metallicity for [Fe/H] = —3, 0, and +1, all for T c g = 
5,000 K and log 5= 5.0. 
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Fig. 3. — Mass loss estimates for metal-poor stars based on Reimers' formula with a variety of 
efficiency parameters rj. Estimation has been made using stellar evolutionary tracks that were used 
to construct the Yale Isochrones 1996 . In the top two panels, the lines with asterisks are from the 
measured mass of the RR Lyrae stars, Mrr, in the globular clusters assuming the age range of 12 
- 16 Gyrs as shown in Table 5. The M RR for Z = 0.0001 and Z = 0.0004 are best reproduced by 
rj < 0.3 and rj s=s 0.3 - 0.5, respectively. 
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Fig. 4. — Same as Figure 3 but for metal-rich stars. The left and right panels are for AY / AZ = 
2 and 3, respectively. Mass loss as a function of age cannot increase indefinitely as rj increases 
because the mass of the red giants and their core mass also depend on age and metallicity. After 
a critical age, mass loss decreases for this reason. 
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Fig. 5. — Lifetime on the MS and RGB, tu (upper panel), the core mass at the tip of the 
RGB, M core , and the mass loss, AM (lower panel), as a function of initial mass M for Z = Zq . 
AM increases rapidly as M decreases because tu increases as M decreases. Below a certain mass, 
M £3 1.0 in case of Z = Z Q , however, AM decreases as M decreases because M core grows too fast 
to allow more mass loss to occur. A different mass loss efficiency 77, therefore, leads to a different 
final core mass. 
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Fig. 6. — Mass distribution functions assumed in this population synthesis. The model with a = 
0.06 Mq reproduces the width of the HB of Galactic globular clusters (Lee et al. 1990). The models 
with a = 0.01 and 0.99 are the two extreme assumptions that represent a pseudo-delta function 
and a pseudo-constant function for mass distribution. 
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This figure is in a separate gif file. 



Fig. 7. — Model CMD for low metallicity as a function of age and metallicity. The synthetic HBs 
are based on the assumption of (77, x, a) = (1.0, 1.35, 0.06). Note that the HB becomes hotter as 
age increases. A larger circle denotes a larger number of stars, where the total initial mass of the 
model galaxy is 10 12 M . It is important to note that these models are based on rj = 1.0. We 
think that this value is adequate only for metal-rich stars and is much higher than the estimated 
value (« 0.3 - 0.5) for metal-poor stars. As a result, the 15 Gyr old Z = 0.004 model with 77 = 
1.0 (bottom right panel) contains too many blue HB stars compared to 47 Tuc. Nevertheless, 77 
has been kept the same for all models of different ages and metallicities in Figures 7 - 12 in order 
to see only the effect of age and metallicity. The brightest PAGB track (in 15 Gyr-old models) is 
the Kiel group 0.598 M model, and the second and third are their 0.565 and 0.546 M & models, 
respectively. The fourth brightest horizontally crossing track is for the stars that do not have a 
sufficiently large mass to experience helium core flash (Sweigart et al. 1974; see text). 
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Fig. 8. — Same as Figure 7, but for metal-rich populations of Z > 0.01 and AY / AZ = 2. Note 
that metal-rich models have many UV-bright (slow blue phase) stars. 
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Fig. 9. — Same as Figure 7, but for metal-rich populations of Z > 0.01 and AY/AZ= 3. The 
UV-bright phase, the SBP, is even more conspicuous than the case of AY / AZ = 2 in Figure 8. 
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Fig. 10. — Model SED for low metallicity as a function of age and metallicity. These models 
are based on the model CMDs shown in Figure 7. A noticeable "UV upturn" appears at a large 
age. Once again, keep in mind that the metal-poor models overproduce UV flux almost certainly 
because rj = 1.0, a factor of 2 - 3 overestimation for metal-poor stars, has been used in this figure. 
The y-axis is in an arbitrary scale. 
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Fig. 11.— Same as Figure 10, but for Z > 0.01 and AY/AZ= 2. The far-UV flux relative to the 
near-UV flux is generally higher than the metal-poor models in Figure 10. 
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Fig. 12. — Same as Figure 11, but for AY/AZ = 3. A strong UV flux develops more quickly when 
AY / AZ is higher because the UV-bright shell helium-burning phase (the SBP) is more significant 
when helium abundance is higher ( Dorman et al. 1993] ; YDK). 
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Fig. 13. — The definitions of central and shell helium-burning stages. The crosses and squares are 
for central and shell helium-burning stages, respectively. All symbols denote time interval of 10 
Myr. 
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Fig. 14. — Light contribution from various evolutionary stages at different ages. The upper and 
the lower panels are for 5 and 15 Gyrs old models, respectively, both with (Z, Y, rj, x, a) = (0.02, 
0.29, 1.0, 1.35, 0.06). At an early age, the most UV light comes from PAGB stars, while the UV 
flux is almost negligible, as seen in Figure 12. The lower panel shows that much of the UV light 
at 15 Gyr, especially in the far-UV, originates from shell helium-burning (evolved HB) stars in the 
UV-bright phase (SBP) discussed in YDK. 
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Fig. 15. — Effects of age and metallicity on the model two-color diagram, m(1500)-V vs m(2500)- 
V. Each line is for one metallicity but covering 1-25 Gyrs of age. From the smallest symbol to 
the largest, symbols denote 1, 5, 10, 15, 20, and 25 Gyr old models. In order to show only the 
effects of age and metallicity, other parameters are fixed, i.e. (rj, a, x) = (0.5, 0.06, 1.35). The 
model colors are defined in the text, and the cluster data are from Table 1 of Dorman et al. (1995). 
Filled circles and filled squares are metal-poor cluster data (Z < 0.002) and metal-rich cluster data 
(Z > 0.002), respectively. Most data points are matched by the models with the age derived from 
the MS turn-off in the CMD (approximately 15 Gyr). Because of the opacity effect, the more 
metal-rich, the redder the UV-to-V colors. The issue about the true rj is discussed in Section 4.2. 
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Fig. 16. — Same as Figure 15, but for metal-rich systems in the m(1500)-y vs m(1500)-m(2500) 
diagram. In order to show only the effects of age and metallicity, other parameters are fixed, i.e. 
(77, a, x) = (0.7, 0.06, 1.35). The galaxy data are from Table 2 of Dorman et al. (1995). If the 
majority of stars in elliptical galaxies are approximately 1-2 Zq , then the r? = 0.7 models suggest 
that elliptical galaxies are about 15 Gyrs old. 
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Fig. 17. — Light contribution from various evolutionary stages for different metallicities. The 
upper and lower panels are for a metal-poor model with (Z, Y) = (0.0004, 0.23), and a metal-rich 
model with (0.04, 0.35), respectively, both with (age, 77, x, a) = (15 Gyr, 1.0, 1.35, 0.06). Most UV 
light comes from the central helium-burning stars (HB stars) in the metal-poor model. But, a larger 
amount of UV light comes from the highly evolved, shell helium-burning stars in the metal-rich 
case (lower panel) due to the UV-bright phase that is more common in the metal-rich stars. 
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Fig. 18. — Effects of age and metallicity on model two-color diagrams. In order to show only the 
effects of age and metallicity, other parameters are fixed as (AY/AZ , r/, a, x) = (3.0, 1.0, 0.06, 
1.35). If r] ~ 0.3 - 0.5 for metal-poor stars as discussed in the text, the metal-poor models with 
T) = 1.0 shown in this diagram are generating too much UV flux at a given age. The issue about 
the true r\ is discussed in Section 4.2. Note that elliptical galaxies are best matched by metal-rich 
models whereas metal-poor models reproduce the globular clusters well as normally expected. 




Fig. 19.— Effects of 77 on model CMDs and SEDs for (Z, Y, age, a, x) = (0.04, 0.31, 15 Gyr, 0.06, 
1.35). F\ is in arbitrary but consistent unit. It is clear that a higher mass loss efficiency leads to a 
larger number of hot stars and a higher UV flux. Various studies favor a high r/, namely, r] > 0.7 for 
Z > Z & . If this is true, old, metal-rich galaxies can exhibit a strong UV upturn within a Hubble 
time. 
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Fig. 20.— Effects of r] on model CMDs and SEDs for (Z, Y, age, a, x) = (0.0004, 0.23, 15 Gyr, 
0.06, 1.35). F\ is in arbitrary but consistent unit. Unless r/ were as high as unity, an unlikely high 
value for metal-poor stars, metal-poor populations would fail to produce a strong UV upturn. 
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Fig. 21. — Effects of 77 on model two-color diagrams, m(1500)-V vs m(2500)-V. See the Figure 15 
caption for model and symbol descriptions. In order to show only the effects of 77, other parameters 
are fixed, i.e. (a, x) = (0.06, 1.35). Data for metal-poor globular clusters (filled circles) and 
relatively metal-rich globular clusters (filled squares) are shown. Note that 77 = 0.3 models, which 
reasonably reproduce the colors of metal-poor clusters, fail to match the colors of the metal-rich 
globular clusters (e.g. NGC 6388, Z « 0.006) at an acceptable age; the r\ = 0.3 model suggests an 
age of about 24 Gyrs for NGC 6388. 
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Fig. 22. — Same as Figure 21, but for metal-rich systems on the m(1500)-V vs m(1500)-m(2500) 
plane. The UV strength of elliptical galaxies cannot be reproduced by low ij models. If elliptical 
galaxies are approximately 15 Gyrs old, this suggests rj ^ 0.7 for metal-rich stars. 





i i i | i i i | i i 


1 




\ This f 


igure is in a separate g 


if file. 













Fig. 23. — An example showing the effects of AY / AZ on model CMDs and SEDs for (Z, age, rj, 
a, x) = (0.02, 13 Gyr, 0.06, 1.0, 1.35). F\ is in arbitrary but consistent unit. A model with a 
higher AY j AZ shows a higher UV flux because (1) HB stars with a higher helium abundance, Y, 
become UV-bright more easily, and (2) stars with a higher Y evolve faster on MS and RGB when 
Z ^ Z Q . Thus, a galaxy with a higher Ay j AZ contains a larger number of low- mass (hot) HB 
stars at a fixed age. 
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Fig. 24. — Effects of AY / AZ on model two-color diagrams. Models are for (Z, 77, a, x) = (0.02, 1.0, 
0.06, 1.35). Higher- AY j AZ models generate UV flux more quickly because stars with higher helium 
abundance evolve faster and because the UV-bright phase of metal-rich HB stars is more significant 
as helium abundance increases (see YDK). But it is difficult to choose the true AY/AZby UV 
color fitting unless the age is known a priori, because, if AY / AZ is within the range of 2 - 3, the 
model sequences of different values of Ay j AZ look very much alike. Symbols are the same as in 
the Figure 15 except that 13 Gyr-old models (filled symbols) are also marked to be compared with 
Figure 23. 
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Fig. 25.— Effects of the IMF slope x on model CMDs and SEDs for (Z, Y, age, a, rj) = (0.02, 
0.29, 15 Gyr, 0.06, 1.0). F\ is in arbitrary but consistent unit. A model SED is not very sensitive 
to the IMF slope x. 
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Fig. 26. — Effects of the IMF slope x on model two-color diagrams. Models are for (Z, Y, rj, a) 
= (0.02, 0.29, 1.0, 0.06). Note that the IMF slope has little effect on the UV-to-V colors. See the 
Figure 15 caption for details about symbols. 
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Fig. 27. — Effects of the gaussian mass dispersion factor a on model CMDs and SEDs for (Z, Y, 
age, a, x) = (0.02, 0.29, 10 Gyr, 0.06, 1.35). F\ is in arbitrary but consistent unit. When a galaxy 
is so young that a larger fraction of its core helium-burning stars is cool, a larger a causes a higher 
UV flux. 
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Fig. 28. — Same as Figure 27, but for a 15 Gyr model galaxy. When a galaxy is old enough that 
a larger fraction of its core helium-burning stars is hot, a larger a causes a lower UV flux. 
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Fig. 29. — Effects of the mass dispersion parameter, a, on model two-color diagrams. All models 
are for (Z, Y, r], x) = (0.02, 0.29, 1.0, 1.35). It is clear that the model UV-to-V colors are sensitive 
to a which determines the width of the HB. Thus, oversimplistic HB treatments are likely to lead 
to wrong conclusions. See the Figure 15 caption for details about symbols. 
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